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I. INTRODUCTION 


The subject grant (NGR 39-009-162) on technical problems related to 
deployment and retrieval of spinning satellites is at mid-term in a one- 
year extension of the original two-year grant. This is the fifth semi- 
annual status report on the project. During the past six months (1 January 
1972 to 30 June 1972) work has progressed in the two major areas outlined 
in the new statement of work. Appendix A of the fourth progress report, 
January 1972. Quantitative analyses have been developed to the point 
where numerical results can be easily obtained in the optimal transfer 
trajectory problem. This effort should be completed within the next few 
months. Analyses of deployment dynamics- and control of the paired- 
satellite concept (AMAPS) are progressing well and results are expected 
by the end of this year. A summary of accomplishments is offered in 
section III. Individual personnel assignments are cited in section II, 
and expected progress for the next grant period is summarized in section IV. 
Detailed descriptions of accomplishments are left for appendices. 

Interaction with NASA and industry is considered essential to avoid 
duplication of effort and to provide helpful guidance. Recently (June 1) 
the principal investigator visited the Astronics Lab. at NASA-MSFC to 
discuss satellite retrieval problems. Just prior to that Mr. Heinz Fornoff 
presented a seminar at Penn State on the subject of manipulators in space. 
Such interactions have been very helpful in directing research tasks. 
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Dissemination of results is a primary function of such research grants. 
In connection with this grant a paper entitled "The Problem of Docking 
with a Passive Orbiting Object Which Possesses Angular Momentum," (See 
Appendix B, Fourth Progress Report, January 1972) will appear this 
Summer in Astronautical Research 1971 , Proceedings of the 22nd l.A.F. 
Congress, published 1972. In addition, the paper entitled, "Attitude 
Dynamics and Control of an Apogee Motor Assembly with Paired Satellites," 
(See Appendix C, Fourth Progress Report, January 1972) has been published 
in the Journal of Spacecraft and Rockets , June 1972, page 410. 

Notification has just recently been received that an abstract 
entitled, "The Problem of Deploying and Retrieving Spinning Payloads in 
Orbit," has been accepted to the Fifth Symposium on Automatic Control in 
Space, sponsored by the International Federation of Automatic Control 
(IFAC) , to be held June 5-9, 1973 in Genoa, Italy. This abstract is 
included here as Appendix A. 


II. PERSONNEL 

The grant budget currently supports three graduate assistants on 
half-time schedules. This project is currently fully staffed with three 
master of science candidates (E. C. Thoms, R. J. Cenker, and A. A. 
Nadkarni) . Two of these students (Thoms and Cenker) should be completing, 
their theses during the latter part of 1972. These will be published as 
Astronautics Research Reports of the Aerospace Engineering Department of 
Penn State, and submitted to appropriate journals for wider dissemination. 

Basic task assignments presented previously have not changed signifi- 
cantly. However, as progress is made interactions and objectives must be 
modified to reflect current and expected situations. Therefore, new 
project flow charts have been constructed showing current areas of study. 
These are presented in Figure 1. 
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Figure 1. Project Flow Chart 

















Figure 1. (continued) 
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III. PROGRESS TO DATE 

The current one-year renewal of this project has permitted a 
significant increase in forthcoming results related to retrieval and 
deployment technology. Specific areas of concentration include optimal 
transfer trajectories for retrieval, attitude reorientation during 
deployment sequences, and automatic controls for retrieval spacecraft. 

Retrieval package transfer trajectories between the orbiter and 
target can be optimized for minimum fuel or time through the use of 
optimal control techniques established in the literature. The problem 
considered here is nonlinear with two-point boundary values and is 
constrained at the terminal point because of the desired spin-axis 
alignment. Such problems require the use of a computer for solution. 

A detailed discussion of progress on this problem is presented as 
Appendix B. 

The initial work reported in the Journal of Spacecraft article, 
"Attitude Dynamics and Control of an Apogee Motor Assembly with Paired 
Satellites," is continuing with the objective of active spin reorientation. 
Progress to date is concerned with developing relationships for tumbling 
dynamics of non-symmetrical satellites to determine limits on nutation 
angles, ratios of moments of inertia, and energy state. Such relationships 
define the critical point in energy level during passive spin axis transfer 
of an AMAPS type device at which final spin axis direction may be controlled 
most easily. A detailed discussion of this problem is presented in 
Appendix C. 

As the retrieval package operational plan developed several controls 
problems arose in connection with attitude maintenance during the capture 
sequence. Such problems are especially important to more automated 
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missions for retrieval, e.g., during synchronous orbit missions'! There- 
fore, work has continued on synthesizing control laws for the various 
situations of interest. Specifically, considerations include the case 
of general target spin vector misalignment relative to the despin ring 
and control of the ring speed during capture. Control system developments 
for these functions are discussed in Appendix D. 

IV. FUTURE TASKS 

Efforts will continue in the areas of optimal transfer trajectories, 
spin reorientation control, and control system synthesis for retrieval 
package functions. Optimal three-dimensional transfer trajectories will 
be generated for various parking positions of the orbiter and target 
spin orientations. Comparisons will be made with results from the 
conjugate-gradient technique and other methods. Recommendations concern- 
ing the best stand-off position for the orbiter will also be included. 

This task should be complete by the end of 1972. 

The work on spin reorientation will continue through computer 
simulations of various control schemes to insure proper final spin vector 
direction. Methods in which moments of inertia are varied and energy is 
changed will be considered. Appropriate mechanisms will be conceived for 
such applications. This task should also be- completed- by the end of 1972. 

Activities related to automatic operations will continue and, 
hopefully, develop into a phase of a new grant involving studies of self- 
contained high-altitude deployment packages, multiple and automated 
retrieval missions, and high-altitude retrieval techniques. A great deal 
of knowledge relevant to orbital retrieval and payload deployment has been 
gained, and each new area of technology has led to further questions 
about retrieval and deployment. 
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APPENDIX A '* 

Abstract for Fifth Symposium on Automatic Control in Space 
(Genoa, Italy, June 5-9, 1973. Sponsored by IFAC) 

The Problem of Deploying and Retrieving Spinning 
Payloads in Orbit* 

Marshall H. Kaplan 

Associate Professor of Aerospace Engineering 
The Pennsylvania State University 
University Park, Pa. 16802 
U.S.A. 

ABSTRACT 

Two of the primary missions for the space shuttle system will be the 
deployment and retrieval of payloads. Many of these will be spinning for 
stability. A large number of old satellites and other debris in orbit 
are spinning due to initial deployment maneuvers or later mishaps. In 
either case, the presence of passive angular momentum requires special 
techniques which will not endanger the shuttle orbiter vehicle. Separate 
spacecraft, operated from the shuttle, may be used to deploy and retrieve 
spinning objects. A conceptual design of such vehicles has been proposed 
and many of the technical problems have been defined.^" 

The paper proposed here will consider many of the "key" problem 
areas related to dynamics and control of deployment and retrieval in both 
low and high orbits. Topics will include: (1) optimal transfer of a 

retrieval package from shuttle to spinning objects such that approach is 
accomplished along the spin axis, (2) control of spin reorientation during 
high altitude deployment, and (3) automatic stationkeeping of the shuttle 

* 

This work is supported by NASA Grant NGR 39-009-162. 
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relative to a retrieval target. The second item represents an extension 

of work on a special technique for deploying pairs of high altitude sat- 

2 

ellites without the use of an orbit-to-orbit tug. During large angle reor- 
ientation of a passive semi-rigid body, there arises an ambiguity in final 
pointing direction of the stable spin axis. The elimination of this 
ambiguity can be accomplished by an active mechanism within the body. 

This will be discussed in detail and quantitative results presented. 

The optimal transfer problem arises in order to minimize fuel or time 
during the retrieval of a spinning object. In the case of a multiple 
retrieval mission the intermediate spacecraft would have limited fuel for 
maneuvers. In other situations time could be the critical factor. The 
requirement that final approach of the retrieval vehicle be along the object 
spin axis represents an end constraint which significantly complicated the 
problem. Automatic stationkeeping, of the. shuttle during these maneuvers 
is closely related to the transfer trajectories since this represents an 
optimum stand-off position for the orbiter. Considerations for selecting 
a combination of parking positions and transfer trajectories will be included. 

References 

1. Kaplan, M. H. et al. , "Dynamics and Control for Orbital Retrieval 
Operations Using the Space Shuttle," NASA-KSC TR-1113, Vol. 1, May 
1971, pp. 175-201. 

2. Kaplan, M. H. and Beck, N. M. , "Attitude Dynamics and Control of an 
Apogee Motor Assembly with Paired Satellites," Journal of Spa cecraft 
and Rockets , June 1972, pp. 410-415. 
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APPENDIX B 

Optimal Retrieval Transfer Trajectories 
(E. C. Thoms) 


I. Introduction 

During a retrieval mission the shuttle orbiter must initially park 
near the target object, preferrably in the same orbit. Then the retrieval 
package is deployed and transferred to the object such that rendezvous is 
along its spin axis. It is very desirable and possibly essential to have 
an optimized transfer trajectory, either to minimize time or fuel (impulse) 
for the maneuver. The optimization technique used must handle a special 
terminal state constraint, spin axis alignment. This represents a diffi- 
cult problem because it leads to a three-dimensional trajectory (out-of- 
the-plane). The transfer represents a solution to a nonlinear boundary 
value problem with a terminal constraint on the approach direction. 

II. Optimization Methods 

An investigation of various optimization techniques has shown that a 
conjugate-gradient method should yield valid results with a reasonable amount 
of computational time. The conjugate-gradient approach Attempts to 
combine the advantages of both first and second order gradient methods while 
eliminating some- disadvantages. As the - number of' iterations increases 
convergence becomes rapid near the optimum solution. In addition, a 
poor initial guess will yield satisfactory solutions, whereas the. first- 
order method would not.^ - 
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In a first-order gradient technique the direction of movement is made 
from a nominal point along the direction of steepest ascent. Most steepest 
ascent methods tend to give an oscillating approach to the local optimum. 

A conjugate-gradient method utilizes gradients in fixing the direction of 
movement, but does not necessarily move in the direction of steepest ascent. 
In this method searches are performed in a certain set of directions (i.e., 
conjugate directions) chosen so that the optimum of a quadratic function is 
found in a finite number of iterations. Quadratic behavior is observed 
near the optimum in all functions which are expandable in a Taylor series. 
However, in the general search problem the coefficient matrix of the 
quadratic term in the objective function is not known and it is not possible 
to compute conjugate directions in advance. 

Various proposals have been made to generate a set of these directions. 

2 3 

Fletcher and Powell have modified a method proposed by Davidon to generate 

4 - 

a set of conjugate directions using the gradient, and Fletcher and Reeves 
have applied a technique suggested by Hestenes and Stiefel^ in connection 
with the solution of linear equations to find a set of conjugate directions. 
Lasdon, Mitter, and Warren^ have extended the Fletcher-Reeves conjugate 
gradient method to function space problems. 

The basic problem associated with these methods is the generation of 
a set of conjugate directions when the form of the quadratic is unknown. 

It is possible, however, to generate sets of conjugate directions based 
on a knowledge of the gradients. The method of Fletcher and Reeves uses 
this idea and appears to be efficient in search and is easily programmed 
for the digital computer. The Fletcher and Reeves method selects the 
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This method is, therefore, guaranteed to locate the optimum of any quadratic 
function in at most n stages since this equation generates conjugate directions. 
For functions which are not quadratic, convergence cannot be guaranteed in 
n steps. Therefore, a test for convergence is needed. The directions of 
movement are related to some quadratic approximation, and depending on the 
type of objective function, convergence to the vicinity of the minimum may 
be fast or slow. However, once near the optimum, the method converges 
rapidly, because most functions tend to behave quadratically in this region. 
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IV. References for Appendix B 

1. Bryson, A. E. and Ho, Y. C. , Applied Optimal Control . Gin and Co., 
Waltham, Mass., 1969. 

2. Fletcher, R. and Powell, M. J. D. , "A Rapidly Convergent Descent Method 
for Minimization," British Computer Journal . June 1963, pp. 163-168. 

3. Davidon, W. C. , "Variable Metric Method for Minimization," A. E. C. 

R and D Rept. , ANL-5990 (Rev.) 1959. 

4. Fletcher, R and Reeves, C. M. , Function Minimization by Conjugate 

. Gradients," British Computer Journal . July 1964, pp. 149-154. 



13 


5. Hestenes, M. R. and Stiefel, E. , "Method of Conjugate Gradients for 
Solving Linear Systems," Rept. 1654, National Bureau of Standards, 1952. 

6. Lasdon, L. S., Mitter, S. K. and Warren, A. D. , "The Conjugate Gradient 
Method for Optimal Control Problems," IEEE Transactions on Automatic 
Control . Vol. AC-12, No. 2, April 1967, pp. 132-138. 
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APPENDIX C 

Deployment Attitude Dynamics and Control 
(R. J. Cenker) 

I. Introduction 

It has been established that a satellite initially rotating about its 
minor axis (axis of minimum moment of inertia) will experience a transfer 
of spin axis due to energy dissipation. The final state is one in which 
stable spin about the major axis (axis of maximum moment of inertia) is 
established. However, an ambiguity in pointing direction occurs. The 
spin vector has two possible orientations with respect to the body axes, 

180° away from each other. Elimination of this ambiguity through active 
control is sometimes essential. Application of the AMAPS technique is a 
particular example of interest here. The situation with ATS-V is another 
case where the ambiguity arose to result in wrong spin. Control of this 
ambiguity is the topic of interest here. 

II. Summary of Past Work 

Effects of energy dissipation on satellite dynamics have been experienced 

1 2345 

in actual flights and discussed in several analytical treatments * * * . 

These efforts have, in general, dealt with discussions of the effect of 

dissipation on the precession rates ; and the total time required for a 

satellite to go from one spin state to another. References 2 and 3 describe 

attitude drift for a symmetrical body and for "n" symmetrical bodies, 

respectively. Reference 4 discusses the motion of a non- symmetric body in 

body fixed and inertial coordinates, while Reference 5 illustrates several 
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methods of modeling the energy dissipation process and describes effects 
on non-symmetric bodies. 

Control of this attitude drift is discussed in References 6,7, and 8. 
Axisymmetric spacecraft control is proposed in Reference 6 , in which a 
satellite can be maintained in a configuration other than spinning about 
the axis of maximum inertia. This is done by addition of energy into the 
system to compensate for energy dissipation. Reference 7 touches briefly 
on non-symmetric satellite precession control with controlled energy dissi- 
pation, and presents numerical examples for symmetric satellites. Control 
of a non-symmetric satellite is covered in Reference 8, but only for small 
perturbation angles, i.e., active energy dissipation is used to return the 
satellite to the desired stable position when it has been perturbed. Thus, 
none of the control systems proposed to date considers the question of 
whether the spin vector will align itself with the positive or negative 
axis of maximum inertia when the spin axis is originally offset by a large 
amount. 

III. Analysis 

Constant energy motion of a non-symmetric, rigid body, in a torque 

free environment is geometrically described by Poinsct's motion. This is 

the motion of an ellipsoid (having axes which are proportional to the 

moments of inertia of the body) that rolls without slipping on an "invariable 

plane" while its center is a fixed distance from the plane (see Reference 5, 
x /2F 

Appendix C). This distance is given by — ~ — when T is the rotational kinetic 

n 

energy and H is the angular momentum. The allowable paths for this motion 
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in the body fixed coordinate frames are closed paths called polhodes , each 

polhode having a fixed discrete energy level. 

A logicial extension of this motion for dissipative bodies is a change 

in the distance from the center of the ellipsoid to the invariable plane, 

thus a continuous change in polhode curves. Figures 2,3, and 4 illustrate 

the above concepts, 5 along with a geometric interpretation of the critical 

point in -the tumbling motion, i.e. , whether w has a positive or negative z 

-2 

component when the energy reaches H / 2I 2 determines the final spin orientation 
Since the nutation angle 0 (position of the z axis with respect to the 
final rotation vector, and therefore the angular momentum vector) is the 
parameter of concern, equations giving 0 as a function of time and/or 
rotational energy are desirable. Although the equations of motion for a 
non-symmetric body are elliptical in nature and, therefore, not readily 
reduced to an analytical solution, limits on the nutation angle can be 
developed as straightforward functions of the energy state. Consider first 
the case with the rotational energy below the separatrix energy (T < H 2 /2I 2 ). 
Using the symbols of Reference 5 the following relationships are developed 

, _ z « 

*' (T.-IO 

. El O x = t - 

where 

6^ = upper limit on nutation angle 

9^ - lower limit on nutation angle 
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With the energy above the separatrix energy (T > H the following 

parallel development can be made. Figure 5 shows the inertia ellipsoid 
with the z axis in one of its extreme positions. As both extreme 
positions are similar in nature (u^ = 0 and y axis out of the plane of 
the paper) the upper and lower limits are 90° +0 , respectively. 
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The above relations take into account the fact that at this point in the motion 
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Bearing in mind the particular energy and momentum conditions at this point, 
the final term in brackets equals zero; therefore. 
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where is the maximum possible energy state for a given angular 

2 

momentum. Thus, T = H /2T, . Making the above substitutions , the 
* max 1 ° \ 

bounds on the nutation angle reduce to the following forms: 

2 

a. Motion physically below the separatrix (T > H /21^) 


J-%- 

/'X ; 3 


o ' / / -~T*- 

^ y /- j[ n 
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b. Motion physically above the separatrix (T < H fll^) 

Q u - - CtW-CA A/W. 

0/ — O/t-CaXw. 

At this stage, two relationships are to be noted: 

1) The intermediate moment of inertia (contained in the term I^) 
has no effect on the nutation angle range below the separatrix, 
and no effect on the lower limit of 0 above the separatrix. 




2) At the separatrix 


T - 




Therefore, the upper limit on 0 is 90°. For a specific body having 
a minimum to maximum inertia ratio of one half ( 1^3 = 0.5) the range of T^ 
is given by 0.5 T^ _< 1.0. Figure 6 shows a plot of the bounds on 0 vs T^, 
assuming that I ^ = * 33 * Therefore, the first pair of equations holds for 
the entire range of T*. Figure 7 shows the bounds on 0 for a body having 
the same * 33 * but with an intermediate ratio, = 0.8. This can be seen 
as the critical point in the motion, as the bounds on 0 divide into two regions: 
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one ending in a final nutation angle of zero (z, to, K pointing in the same 
direction), the other ending in a final 0 of 180° (z. pointing in the opposite 
direction as to and fi). 

Noting the continuity of the minimum values for 9 at the transition 
point, the following relations can be written: 

% * ,0 ‘- 
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This coincides with the earlier observation that I ^ did not affect the 

bounds 0 9 , 0". , and 0. , as the value for I.» has no effect on these 
T> ^b a 1 

portions of the curve. 
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From the above description of tumbling motion it is apparent that 
the critical point in the motion is where the energy satisfies 

T = H 2 /2I 2 


At this point either the nutation angle or the z component of to identifies 
which direction the final spin vector will assume. Final control over 
this direction will be attempted by assuming some realistic energy dissipa- 

.9 

tion rate and either changing the moment of inertia by moving some mass 
within the spacecraft, or altering the energy state by starting an electric 
motor within the satellite. At this point, any changes in the moments of 
inertia are expected to use symmetrical motions of masses to avoid shifting 
of the center of mass of the satellites, simplifying the resulting equations 
of motion. Present plans also call for lining up any control motor with the 
z axis, as this component of the rotation appears the most critical. 


IV. Conclusions 

Results to this point yield a critical state at the point where 
2 

T = H 17-^2’ At this point, either the nutation angle or the z component 
of the rotation vector can be used to predict the final spin orientation, 
according to the following criteria. 
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Prior to this point (T = H {21^) the z axis oscillates about 90° from 
the angular momentum vector. After passing this point, it oscillates 
through a steadily decreasing range to either 180° or 0° from fi. Control 
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over this final orientation will be sought rising either available hardware 
on the spacecraft (extendable booms, electric motors) or equipment 
designed especially for control purposes. A computer simulation will 
be developed to demonstrate the feasibility of the scheme, and areas of 
potential improvement or difficulty will be identified. 
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APPENDIX D 

Automatic Control Systems for the Retrieval Package 
(A. A. Nadkami) 

I. Dynamics and Control during Capture with Spin Axis Misalignment 

The analysis used to predict disturbed motion of a spinning object 
due to accidental application of a torque by the ring'*' is extended to 
a more general case of spin-vector misalignment. Increased divergence 
of the composite body (target plus ring together) during despin is 
considered. Associate control systems and limitations and misalignment 
are being studied in addition to an appropriate damper for stabilizing 
the dynamics. 

At least a small alignment error can always .be ...expected during the 

retrieval sequence. The general situation considered is illustrated in 

Figure 8. It is assumed that U) g - u) r , i.e., the two spin vectors are 

approximately synchronized before capture. To calculate the resultant 

2 

angular momentum of the total system after capture use 

- • _ « 

W * Ho W, -h X rrA 0 -V rq y An, rq (1) 

where H q and represent the angular momenta of the reference body 
and satellite about their respective centers of mass. 

"Referring to Figure 8, the total angular momentum about the composite 
center of mass C is 

4+ t - H 0 + * p X /*, X r^\ fi (2) 
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The total angular velocity of the composite body (satellite and 
despin package with ring) is given by 
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Now, using generalized momenta 
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we obtain Lagrange's equations of motion 
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For free motion 
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By substitution of (17) into (14), it is seen that sin 0 is a Jacobian 
elliptic function of a linear function of time. Comparison of the 
solutions of this type yields the modulus K of the elliptic functions. 
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Solving for the Eulerian angles 9, 4* > and 4> in terms of elliptic functions 
we have. 
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(Lo^QA'^y 1 -— iA=. — f_ c (21) 

(c^ V- + 7 > *4-_ (\ -z ^ u«- 2 ^;(r -t z +. . ^ 


6r - 


(Ol^X y- c^^uLL ^ 


(22) 



The constant of integration may be calculated from the initial conditions. 
Thus , the motion of the composite body is given in terms of the Eulerian 
angles as functions of time. 

The angular speed of the composite body will be maintained provided the 
spin rate of the design ring is constant. After capture of the satellite. 



the retrieval package control system switches off the despin motor. 

The ring and satellite are then despun to zero angular velocity before 
the retrieval package returns to the shuttle. Assuming that the despinning 
operation takes a finite amount of time, we can assume a linear decay 


law for the non-constant spin rate., 

1o 


+ ' 


(li- CLjt') 


(24) 


From Euler's equations of motion, assuming that the body is approxi- 

4 

mately mass symmetric about the rotational axis, we may write, 

I 

- I (i + -^r)+ (25) 


and 




f (26) 

where use has been made of cos 0-1 and ip sin 0 << p. Combining equations 
(25) and (26), yields 






(27) 


Equation (27) governs the pitching and yawing motions, with respect 
to the space axis, of the rotationally symmetric bodies with non-constant 
spin rates. 

The general form of the solution to (27) is 




= e 



+ + 1 0 (28) 

Equations for 0 and ip could be obtained by separation of the real and 

4 

the imaginary parts of equation (28) 
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II. Despin Ring Speed Control 

Once the retrieval package is aligned with the spin axis of a target 
satellite the despin ring is spun up to the same angular rate as the target. 
After the speeds are synchronized capture arms extend inward to the 
satellite. As these arms move inward the inertia characteristics of 
the ring are changed and its speed will increase. A control system is 
desired to maintain the speed as constant as possible while the inertia 
changes. An initial system design with feedback proved to be unstable 
with the values of the variables chosen as being typical. A compensated 
system was then design with another feedback loop which greatly increased 
stability. 

The voltage applied to the field is obtained from the output 
of an amplifier in low-power operation or from a dc generator when 
greater power is required. In the field circuit, the resistance of the 
windings is R^, and the inductance is designated by L^. 

The torque T developed by a motor is proportional to the product 
of the armature current I & and the magnetic flux cj> of the field: 

T- K. 

where K q is a constant for any: motor and depends upon the total number 
of armature conductors , the number of poles , etc. 

A typical curve of flux 4> versus field current 1^ is shown in 
Figure 9. When the field current becomes great enough to cause 
the iron to saturate, the flux <p no longer increases linearly with 
the current. Motors used in control systems usually., operated over 
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linear portion of this curve, in which case: 

<jf> ■= Ki.^4 

where is the slope of the linear portion of the curve shown in Figure 
9. Therefore, we may write: 

T- 

where 

tC, = KK*. . 

If the moment of inertia of the armature is J, the coefficient 
of viscous friction B v , and the load torque T^, then from a summation 
of torques acting on the armature it follows that 

where 1 

^- 7 i. 


The equation for the field current 1^ is obtained from the equivalent 
field circuit 

C-f 






where '~Ec r ~ — is the time constant of the field circuit. 
V *Sr 

The circuit equation for the armature is then 


(e„-Kc&DKi - R*0 + 

- 

' K4C 1 +T +V S )I + 



38 . 


where 




and 





(command input) 

voltage signal from tachometer 


The transfer function is defined as follows where 






Kj ~~ \ __ — p ' 

I + T^>> J L 3 



&c. _ KT KV Ka* 

+ k; KiK m t 


In Laplace form with all initial conditions equal to zero we have 


<kcp 
4 CO 


KiKrK»t/ T*.? 


5V( 


13 




s v f IzAs + it.kik^C 

^st+' -ZfT 


For the disturbance T , the transfer function has the form 

J-i 

__ 4(M 

TuCO sO 
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The characteristic equationsof the above transfer functions rendered 
the system to be marginally stable at best for all values of the varied 
parameters. 

Compensation was added to the system and the resulting block diagram 
appears as in Figure 10. 

The characteristic equation for this sytem with 'Vj J . = ^z. and 

negligable is : 

sV -b s' 3 4- = o 

where 

|— I K^I Kt-K/y^ C 

___ 

K) VCr K™ G- 

T- 

The transfer functions thus become 

4.^ 

and 

4 ££> 




4 r,. 


K 


G , 




Figure 10. Compensated System 
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The block diagram for the compensated system appears in Figure 10. 
Using the characteristic equation 


5 




S 3 + R 


ti 


SV G 


't, 


+. 

2 V 


G_ 


root-locus diagrams were dran for what was felt to be various typical 
values of the parameters T 2 and H. 

A good indication of the effect of the ratio °^ ta ^ ne< ^ by 

constructing the root-locus plot for a small ratio such as ^ 1/^2 = 
and for a large ratio such as = 10.0. 

Tables 1 and 2 show the scheme used for the root- locus plots. 

The most satisfactory plot is shown in Figure 11. 

Here we have a fast T 2 and a relatively large range on the value 
of G for stability. The loci were calculated for values of G up to 
200 and the system is completely stable. A value of G = 100.0 gives 
a damping ratio on thb inner loop of ^ = 0.95 with 3 = 18°. The natural 
frequency is = 5.2631. 

The time response of the system appears as 




S' 7. t e. (j • (a Z- -+- 1 s 1 , $ eN ) 

+ (&.v 4- ir5.S" QV ) 


The disturbance time response is 


& c C£) 


-.437 

“S' 
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-s 


(/-Ctb Jtl.l*/} 

+ i3s.r) 


TlC*) 
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X h =10.0 
H 



10 

20 

40 

0.05 

U 

S 

S 

0.10 

S 

S 

S 

0.50 

s 

S 

S 

1.00 

u 

U 

U 


U = unstable 
S = satisfactory 

TABLE 2. Root-Locus Scheme for Large T ^/ T 2 
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The steady state value for a step reference input are (G ^ ^ 

and for a step disturbance input C^ss- -'/o oo>T . 

The disturbance torque T is due to the changing inertia char- 

Li 

acteristics of the despin ring. Once the ring is spun up to the 
speed of the target satellite the capture arms move inward and decrease 
the inertia of the ring. The angular speed of the ring will increase 


and .be out -of sync, -with that of the target . The inertia of the ring 

2 

with the arms fully retracted is 15160.5 lb ffl -ft and with the arms fully 

2 

extended 10360.5 lb -ft . One might assume the torque T to be a 

in Xj 

linear function of distance. This may be approximated by a ramp function. 
But at a certain point the arms will stop and the ramp function will 
level off to a constant value. 

In order to resolve this we may regard the input as being the sum 

of separate functions as is illustrated in Figure 12. The sum of the 

ramp function which begins at t = 0 and the equal but opposite ramp 

function which begins at t = t Q is seen to yield the function of Figure 

2 

12. The transform of the first ramp function of Figure 12 is x/S , and 

2 -t s 

that for the delayed ramp function is - (x/S )e o . Thus the transform 

T (S) for the input is 
L 



Now, 

Q c O'} * 

(V> |0Srt7.(.>f)(s z + US+lt.-T) V 1 ' 

(_SV /OS 4- Z7.l^ (_s'-+ l«S 
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•f* 


* 


— t s 

Because of the delaying factor e o the second term on the right- 


hand side of the above equation should be ignored for t < t Q . 
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